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Q\ \ Abstract 

We use spin-coherent states as a time-dependent variational ansatz for a semiclassical 
description of a large family of Heisenberg models. In addition to common approaches 
we also evaluate the square variance (H 2 ) — (H) 2 of the Hamiltonian in terms of coherent 
states. This quantity turns out to have a natural interpretation with respect to time- 
dependent solutions of the equations of motion and allows for an estimate of quantum 
fluctuations in a semiclassical regime. The general results are applied to solitons, instan- 
tons and vortices in several one- and two-dimensional models. 
PACS numbers: 75.10.Hk, 75.10.Jm, 03.65. Sq 
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1 Introduction 



There has been a lot of study on classical spin systems in the last decades, commonly using 
a continuum description, for a review see Kosevich et al. 1990. Some one-dimensional mod- 
els have been identified as integrable and treated by means of the inverse scattering method 
(Takhtajan 1977, Fogedby 1980, Sklyanin 1979). A common feature of such nonlinear evolu- 
tion equations is the existence of localized solutions like solitary waves and, in the integrable 
case, mathematical solitons. 
In quantum spin systems analytical treatments appear to be restricted to one-dimensional 
models. For the spin-^-Heisenberg chain of N lattice sites Bethe (1931) managed to reduce 
the diagonalisation of the Hamiltonian acting on a 2 Ar -dimensional Hilbert space to the prob- 
lem of solving N coupled nonlinear algebraic equations, which are well understood in the 
thermodynamic limit. By this Bethe ansatz and its algebraic formulation many properties 
of the model such as thermodynamic quantities could be extracted. Nevertheless, the eigen- 
states of the Hamiltonian themselves are obtained only rather implicitly and are naturally 
translationally symmetric, i .e. non localized. Thus, a quantum analogon to a classical time- 
dependent soliton, which should not be an eigenstate of energy but a superposition of them, 
has not been found yet. 

So it is desirable to give a description of solitons in the semiclassical regime, i.e. for large 
but finite values of the spin length. Considerable work in this direction has been done using 
different kinds of bosonization and tensor products of coherent oscillator states as variational 
ansatz, see e. g. De Azevedo et al. 1982, Ferrer and Pozo 1988, Kapor et al. 1990. Disadvan- 
tages of these methods are difficulties in treating infinite series of nontrivial operator products 
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and, more fundamental, the problem of mapping the finite dimensional Hilbert space of a spin 
onto an infinite dimensional bosonic space. 

An alternative to such approaches are spin-coherent states as introduced by Radcliffe (1971). 
Balakrishnan, Bishop (1985, 89) and Balakrishnan et al. (1990) used these objects for treat- 
ing ferromagnetic Heisenberg chains and evaluated quantum corrections to classical soliton 
dispersion laws. In detail, their way of proceeding has been subject to some criticism by 
Haldane (1986), but the general findings concerning soliton stability with respect to quantum 
fluctuation seem to be valid and will be confirmed in the present work. Moreover, Frahm and 
Holyst (1989) suggested an extension of spin-coherent states introducing a squeezing similar 
to bosonic states. Unfortunately, this approach leads for generic models to very complicated 
equations of motion for the classical angular variables and the additional squeeze parameter, 
that cannot be solved exactly for nontrivial cases. 

In this paper we follow the above authors and use (unsqueezed) spin-coherent states as varia- 
tional ansatz for a large family of Heisenberg models. In addition to the quantities considered 
in the above references we also calculate the quantum mechanical variance of the energy. This 
quantity can either be used as a test of the validity of the approach; on the other hand, the 
general expression obtained has a natural physical interpretation concerning energy fluctua- 
tions in time-dependent spin structures. Finally we apply our results to several prominent 
solitary solutions of spin models. 



2 Coherent states 

In the Hilbert space of a spin of length S we define a spin-coherent state {S;^, </?) by the 
equation 

S0 )V> ■ S\ S; #,</>) = HS \S; r d,(p) (1) 

for the direction sqw = (sin i9 cos (/?, sin lysine/?, cos$). In the usual basis of eigenstates of S z 
(S z | to) = Tim \ m)) these states can be expressed as 

\Si*M = m?) \S) = £ ( 2 f) 2 (cos (^y^ (sin (^)) n e^ n -V \S - n) (2) 

with 



n v n J V \2 

?1=0 \ / 



Ufa if) = exp exp • (3) 

Clearly, we have (S; <p\S \S;$, ip) = fiSs$ !lfi , but in the expectation values of higher products 
of spin components contributions of different order in the spin length 5* arise. A list of diagonal 
elements useful for the following is given appendix [A], where we also demonstrate the classical 
limit of such expressions. 

The spin-coherent states have the minimum uncertaincy product 

A(ei • §) A(e 2 • §) = ^A(e 3 • £) (4) 

with A(O) being the variance of an operator O and ei,e*2,e3 an orthonormal system. Further- 
more, they fullfill the (over-) completeness relation 

25 + 1 f 2 * 



4tt 



dtp d$sm$\S;$,ip){S;$,<p\ = 1 (5) 
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Nevertheless, it should be mentioned that for an arbitrary linear combination of coherent states 
a direction solving eq. (||) cannot be found. Therefore the spin-coherent states do not 
form a subspace of the Hilbert space, but a submanifold diffeomorphic to the two-dimensional 
unit sphere. Only in the case S = | every (normalized) state vector can be identified as a 
coherent state. 



3 Time— dependent variational method 

Let us consider the following class of one-dimensional spin models 



n 



^ ^ [o^n^n+l + (5iS l n S l n + JiS^ 
n i=l 



(6) 



The index i runs over the three spatial directions x, y, z. Each lattice site labelled by n carries a 
spin of uniform length S. The Hamiltonian includes an anisotropic exchange coupling between 
nearest neighbors, a local anisotropy and a magnetic field. The parameters on, ft, 7$ may be 
chosen arbitrarily, in particular with respect to their sign. 
We now employ the direct product of coherent states 



\m) 



\S;$ n ,<p n ) 



(7) 



as a time-dependent variational ansatz, i.e. we assume the time evolution of the state \tp(t)) 
under the above Hamiltonian to be given in terms of time-dependent functions # n (i), (p n (t)- 
With (•) denoting an expectation value within (Q), s n = (sini? n cos <p n , sm$ n sin ip n , cos # n ) 
and cartesian directions e l , we have for TC 



(H) = (hS) 2 £ fl [«« fa • e l ) (4+1 • e l ) + (l - -L) ft (a n • + g (s n ■ e 



n i=l 



(8) 



and from the Heisenberg equations of motion ^.S n = — ^[S n ,H] : 



—s n = {hS)Y, ot t (s n x e 4 ) ((s„-i + s n +i) • e*) + 2 ( 1 - ^ J ft (s„ x e*) (s n • e*) 

i=l 



0) 



The quantity becomes the classical spin length in the limit S — > oo, fi — > 0, hS = const. 
Eqs. (||), @ are identical to the classical energy and the Landau-Lifshitz-equation up to a 
renormalisation of the local anisotropy by a factor (1 — reflecting the fact that this term 
does not contribute to the dynamics for S = 5. Therefore the variational ansatz provides 
semiclassical corrections to the equation of motion and reproduces the correct classical limit. 
This property of spin-coherent states has been found by several authors before with respect 
to particular spin models, see e. g. the references given in the introduction, and is proved in 
appendix |B| for an arbitrary Hamiltonian. 

Next we examine the square variance of the Hamiltonian, i. e. (H 2 ) — (H) 2 - This quantity 
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is nonzero only in the quantum case and (as well as (H)) strictly an invariant of the system, 
whatever the exact quantum mechanical time evolution of the state ([?]) will be. After extensive 
algebra one ends up with the following expression: 

(H 2 ) - (H) 2 = 

(^) 4 E f 



\i=l 

+2 (l - -L) 0i (s n x gi) ( 4 • r) + ^ (4 x e*) ] ) 
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uj'=l 



((V„-i x e J ) • (V n _i x e J ) + (s n +i x e 4 ) • (s n+1 x e J )) 
- (s n • x e j )) ((s n -i + • x e j ))) 

+ 2 E [ f 1 " 5s) ((*« x e ") 2 X " (*» " ( e " X ^)) 



(10) 



The square variance of the energy consists essentially of two contributions being of order ^ 
and -gj respectively. In the above summation over the lattice sites the squared expression in 
the term of order ^ can be recognized as the r.h.s. of eq. @. Thus, we have 

Within our variational approach the leading order of the quantum fluctuations of the energy 
is purely due to the time-dependence of the state vector. On the other hand, for a quantum 
state which has a nontrivial time evolution and is consequently not an eigenstate of the Hamil- 
tonian, the energy must definetely have a finite uncertaincy. Following this observation the 
leading term in (11) is certainly not an artifact of the ansatz (0), but it is a physically relevant 



expression for the energy fluctuations in time-dependent semiclassical spin structures. So we 
have found good evidence that the ansatz of spin-coherent states does not only reproduce 
the classical limit but is still meaningfull for a semiclassical description. The contributions of 
order in ([R]) cannot be interpretated in a general way and should be studied in the context 
of particular models. We will see that these terms can often be considered as a criterion for 
the validity of the variational method. 

For brevity we have concentrated in (^) on an one-dimensional model. For higher dimensions 
one simply has to infer the appropriate number of neighbors of each lattice site in the summa- 
tions. The result can always be written in the form (|Tl|). Moreover, it is also straightforeward 



to see that the result (11) is still valid in the case of exchange couplings of longer range. 
Thus, our above findings apply to a large variety of (ferromagnetic or antiferromagnetic) spin 
models in arbitrary dimension. This will be illustrated in the next section with respect to 
several one- and two-dimensional ferromagnetic models. 
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4 Application to solitons, instantons and vortices 



We now calculate the quantum mechanical energy uncertaincy for solitary solutions in differ- 
ent Heisenberg models. We denote the leading order in ^ by f^i (cf. eq. (|ll|)), the remaining 
contributions by U 2 , i.e. (H 2 ) - (H) 2 = fii + Vt 2 - 

(i) A ferromagnetic Heisenberg chain with isotropic exchange coupling and an external field 
is given by 



n 

with J > 0. Eqs. read: 

(Hi) = -(ks) 2 jy: 



B 

ns 



4rS n = (TiS)J [ S n X (s n -i + S n+ i) +s n X ^ ] 



dt 



ns 



(12) 

(13) 
(14) 



i 



^ \ 2 



B 



7777 { S n X (s n -l + S n+ l) + S n X — 



2S 



ns 



g^2 ~~ SnSn-l) 2 + (1 ~ S*n4+l) 2 ) 



(15) 



We choose appropriate units with J = 1, foS 1 = 1 and the lattice spacing a = 1 and calculate 
the expectation value of (12) in the usual continuum approximation for an infinite system 



(16) 



where we have put the magnetic field in z direction, substracted the ground state energy and 
introduced the canonical conjugate fields p = cost?, q = <p; £ denotes the spatial variable in 
chain direction. The equations of motion (|l4]) read in continuum approximation: 



d t q 



1 



p 



1 — p 2 * (1 — p 2 ) 2 

„2\ q2 



(P&) 2 -p{diqf-B 



d t p = (l-p)d^q- 2p(d^p)d i q 
These equations have well-known soliton solutions (Tjon and Wright 1977): 



with the soliton width 



(B + u)T 2 



sech 



2 ( t ~ vt ~ Co 



qo + cut + ^- (£, — vt — £o) + tan 1 ( — tanh 



C-vt-Cc 



r 



i 



> o 



B + oj 



(17) 
(18) 

(19) 
(20) 

(21) 
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and £; > Qo being constants. The above solution is a pulse soliton with velocity v and an internal 
frequency oj constrained by Its energy is calculated from (|l~6|): 

4 B 4 

E = - + (22 

The square variance (Wf) — {T~Li) 2 = Q\ + ^2 can be obtained from the continuum version of 
©: 

l/l/, 4m 2 \ 1 / 8u, 2 8w \\ 

n. = s(r( 8 ( 2B + ^ + lT^j-F( 16 + I(irW + ir^jj- (23) 



where (21) has been used. As explained in the previous section the quantity f2i can be 
interpreted as a natural property of a time-dependent spin state in the semiclassical regime. 
As a criterion for the importance of quantum effects absent from our ansatz (0) we compare 
f^2 with the soliton energy, i. e. 

Y " (25) 



E S^Jf2V?>{2B + 



We see that higher order quantum corrections to our variational approach are neglegible for 
solitons with large width and consequently small energy. If this condition is not fullfilled the 
spin length S becomes significant, in agreement with the work of Balakrishnan and Bishop 
mentioned in the introduction. 

The continuum approximation is a good description of the system if the classical fields p, q 
vary only weakly on the typical length scale. This means that the soliton width should be 
significantly larger than the lattice spacing a = 1, e.g. T > 10. Thus, continuum and semi- 
classical approximation work well in a consistent area of the soliton parameters. 



(u) Next we consider a ferromagnetic Heisenberg chain with a biaxial local anisotropy: 



n 2 = -jJ2 



SnS n +l + T~xS n S n + T z S n S n 



(26) 



The square variance of the energy is given by 
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- (H 2 ) 2 



i 

25 



s n x (s n _i + S n+ l) 

)2 

+ 8^2 ^ ~~ SnSn-lf + (1 ~ Sn^n+l) 2 

+2 (l - ^) (r 2 (4 x e*) 4 + r| (s n x e 2 ) 4 

+2tvt 2 ((4 • e x ) (s n • e z ) - (4 • (e x x e 2 )) 2 )) 



(27) 



Proceeding as above one can derive the following solutions of the continuum model (see e. g. 
Kosevich et al. 1990): 

p(e,t)=sign(t;)tanh ^~^ + g0 ^ , q&t) = q (28) 

with 

r = I v 2 = (rxsin(2g )) 2 (l- ^) 

^2(1-^)^-^0082 g ) ' 2(r z -TVCOS 2 g ) 

£o, <?o being constants and t z — t x cos 2 go > 0. The eqs. ( p8| ) describe moving kink solitons 
parametrized by go with energy E = S above the ground state. For finite tv the velocity u 
vanishes for g £ {0, tt} and takes its maximum i> 2 nax = |tv cos(2g*)| at go = Q* with 



cos(2g*) = 2— - 1 - sign(— )W(2— - l) 2 - 1 (30) 

TV TV V TV 



Inserting eqs. (28) into the continuum version of ( p7| ) leads to a divergence in arising from 
the term proportional to t 2 . This contribution is almost uniform over the whole system except 
for an area around the center of the kink. Moreover, the same divergence occurs if we evaluate 
0,2 for the classical ground state of the model (p6[). The explanation for this behaviour is 
that, different from the previous model, the quantum mechanical ground state of (^) is not 
described exactly by our ansatz (|7|). Thus, the divergent contribution in order is clearly 
an artifact of our variational approach due to a more complicated structure of the quantum 
mechanical ground state, which does not correspond trivially to the classical ground state 
solution. Note that a similar infinite term in arises for the classical ground state of an 
antiferromagnetic Heisenberg model, e. g. (|12^) with J < 0. Here the divergence in order 
is due to the contribution from the isotropic exchange in the Hamiltonian. 
Since we are interested in physically valid quantum fluctuations in time-dependent excitations 
within a semiclassical approach, we renormalize (Tt 2 ) — (H2} 2 by neglecting 0,2- The remaining 
contribution is calculated easily giving 

3 

Id 2 1 (r x sin(2g )) 2 (l - ^) 2 

Qx = -- = * (- . (31) 

^ 1 * (2(t z -r x cos 2 g )) 2 
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This leading order of energy variance vanishes of course in the static case qo E {0, ^} and 
takes its maximum at go = Q** where 



cos(2<r ) = 22—1- 2sign(^)W(2^ - 1)2 - - (32) 



Obviously, q* and q** do not coincide. 

In the case t x = the kink solitons become static, and the exact quantum ground state is 
again the trivial ferromagnetic state. Thus, Sli = and in no unphysical divergency arises: 



S 2 3V2^ \ 2SJ V V 25 

Comparing y/Tl^ with the energy of the kink excitation, we again conclude that the semiclas- 
sical description should be valid for solitons with large width and low energy. 

{in) The isotropic Heisenberg ferromagnet on a two-dimensional infinite square lattice is 
given by 

(34) 

m,n * ' 

In appropriate units we find in the continuum approximation 

m = U**(^ + (i-?)w 2 ) (35) 



2 J \l-p 2 

Eq. (^) defines a two-dimensional conformal invariant field theory that has been investigated 
by Belavin and Polyakov (1975). The fields p, q parametrize points on the unit sphere S 2 . 
Compactifying the two-dimensional plane by an infinite point we can consider the solutions of 
the Hamiltonian ( |35|) as mappings S 2 — > S 2 . These mappings can be classified in homotopical 
classes characterized by the degree of mapping d with 



If dsP ds^ 1 f 

\ d \ = j d x£a ^ £ ^ uSa d^Td^~ = 4^ J <Mx)d(cos #(£)), ( 36 ) 



which is the number of times that the sphere is covered in the course of mapping. As shown 
by the above authors it holds 

,„ , , , f , 9 1 / dw dw dw dw \ , . 

with w(z,z) = cot -^e ltp and z = x\ + ix^- The minimum in ( |37| ) is realized by instanton 
solutions that are given by arbitrary meromorphic functions w{z) : 

with mi,rii > and ctj 7^ bj for all i,j and the degree given by d = max{J2i m i,J2i n i}- The 
real scale parameter R does not influence the energy of these static excitations due to the 
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confer mal invar iance of the model. 

Let us now consider a single instanton simply given by w = (Jj)" , d= \ or 

\z\ 2 » - R 2 » 



P 



\z\^ + R 2 ' 

Clearly we have fli = 0, and for Q2 one obtains 

1 2vr 2 v 2 - 1 



- , q = v arg(z) 



S 2 R 2 3 sin^ 



\v\ > 1, 



(39) 



(40) 



and = S i R i 4f fo r M = 1- I n contrast to the energy of the continuum model this quantity 
is not independent of R, but scales with Again, for the semiclassical description to be 
valid, y/Tfo should be small compared with the energy (H3) = 47rd. This can be achieved for 
arbitrarily high energies provided that the width of the instanton is sufficiently large. 
Next we consider a solution consisting of two instantons separated by distance 21, 



w(z) 



z-l\» fz + l 



and find 



with 



F(n,v;a) 



2S 2 I 



d 2 z 



R 



I dw |2 
I dz I 



R 



(1 + \w 



2^2 



d 2 z ( \ z ~ a \ \ z + a l ((A* + u ) z + (/" - u ) a ) 



1 + \z — ap^lz + a| 



2/y 



(41) 



(42) 



(43) 



Again, the square variance of the energy scales with jp. Unfortunately, the quantity F(fi, u; a) 
cannot be expressed by elementary functions; e.g. in the case fx = v one obtains after some 
algebra: 

|x| 4/1_4 |x + e| 



F{^ l i-a)=a'^- 2 
with e being an arbitrary unit vector. 



\d 2 x^- 
J (a 



{a-** + \x\ 2 ^) A 

(iv) Finally we consider an anisotropic two-dimensional Heisenberg Hamiltonian: 



(44) 



m,n 



S m ,n I Sm+l,n + <Sm,n+l ) ^^m,n^; 



(45) 



Choosing A > and the lattice lying in the x-y-plane the classical ground state is given by a 
parallel spin configuration in this easy plane. Further we have: 



(Hi) - {H A ) 
{hSfj 2 £ 



)2 
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( ^m,nSm— l,n) "I - (1 ^m,n^m+l,7i) "I - ^m,n^m,n— l) "I - (1 •Sm,n , Sm,n+l) 



+2(1-^) A 2 (4,»xe z ) 4 



(46) 



For sufficiently large A planar vortices are stable excitations (Gouvea et al. 1989). These 
objects are given in a continuum approximation by 

p(x,t) = , q(x,t) = znan -1 f^j (47) 

with energy E = 7rz/ 2 ln(-^), where we have chosen appropriate units as before. The integer 
v is called the vorticity, L is the size of the system and a = 1 is the lattice spacing, which is 
a lower cutoff for the integration. The energy diverges logarithmically with the system size, 
due to the fact that the spin configuration is not parallel far away from the vortex center. In 
this sense, the planar vortex is not localized. 

Inserting such a static planar solution into ( f46| ) obviously gives = and in 0,2 a contribution 
proportional to A 2 that diverges quadratically with growing L. This effect is completely 



analogous to the divergence in the model (26) and has its reason in an unexact description of 
the quantum ground state within our variational ansatz. So we renormalize VL2 by neglecting 
this artificial term, which is identical for any planar solution and therefore in particular not 
sensitive for the vortex (|47]). The remaining expression for the quantum mechanical variance 
of energy is: 

i [w 2 r. r 



In a classical easy-plane model of the above type a topological phase transition, or Kosterlitz- 
Thouless-transition, occurs at a certain temperature Tkt (Kosterlitz and Thouless 1973, 
Kosterlitz 1974). Below this temperature vortices and antivortices (differing by the sign of 
the topological charge v) are bound in pairs, while they become mobile above Tkt- 
Now we use our result ( |48l ) to estimate quantum corrections to the value of the critical tem- 
perature Tkt- We assume that only vortices with \v\ = 1 are present in our system and follow 
some rough approximations in the given references. 

In the free energy F = E — TS the quantum fluctuations will give an additional contribution 
to the entropy S. For a single vortex the classical entropy is S = ln(L 2 ), since the vortex 
can be placed anywhere in the system. In our semiclassical description we have to take into 
account the uncertaincy of the energy: 

»-Ki + f))-." + ¥ + °&) 

The entropy term will dominate in the free energy above a temperature 

which reduces to the classical estimate Tkt = f for L — ► oo. Thus, we conclude that in an 
infinite system there is no quantum correction to the Kosterlitz-Thouless-temperature in first 
order in 4, but a logarithmic finite-size contribution arises in this order. 
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5 Conclusions 



In this work we investigate a large family of spin models in the semiclassical regime with 
respect to quantum fluctuations. Our main result is given by eqs. (|lOl), (|IT|) which provide a 
physically relevant expression of quantum fluctuations of the energy in low orders in ^. In 
particular, the contribution in first order is purely due to the time-dependence of the spin 
configuration in semiclssical description. For a static spin configuration the square variance of 
the energy is of order 4^. The term of first order is a natural property of a state vector with 
a non-trivial time- dependence, which is consequently not an eigenstate of the Hamiltonian. 
The terms of higher order can often be used as criterion of the validity of the semiclassical 
approach. 

These findings are valid for a large variety of spin models in arbitrary dimension. In the 
previous section we have illustrated this by several prominent spin models. 

Acknowledgement: The authors are grateful to Holger Frahm for useful discussions. 
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A Expectation values within coherent states 



Let V[S + , S , S z ] be an arbitrary product of operator-valued spin components. With U given 
in eq. (Bt) we have 



(S;i?,¥>|7>[5 + ,S _ ,S 2 ]|S;i?,¥>) = (S\V[U + S + U, U + S~U, U + S Z U]\S) 



It holds: 



U+S+U = e 1 * cos^ - )S + - sm z - ) S~ + sintfS 



U + S Z U 



- - sin $S + - - sin $S~ + cos 
2 2 



(51) 

(52) 
(53) 



Applying S + , S~ on states typed \S — m) gives a contribution of leading order y/S, while S z 
gives a factor of leading order S. By such arguments it is easily seen that the leading term in 
eq. is proportional to the product of the classical spin components given in s# i¥ ,: 



Jim (S; 0, V \V[S + , S-, S Z ]\S; 0, <p) = P[S d s+ v , S d s^, S d s%J , 



RS = S c! =const. 



(54) 



(S;$,ip\S z S z 


\S;d,<p) = 


(S;$,<p\S + S- + S~S+ 


\S;#M = 


(S;$,<p\S + S z + S Z S + 


\S)#,<p) = 


(S;#,ip\S z S z S z 


\S)#,<p) = 


(S;$,cp\S + S z S + 


\S)#,<p) = 


(S;$,cp\S z S + S z 


\S)#,<p) = 



UiS)- ( cos 2 d + sin 2 •& 



where (hS) = S d = const, is the classical spin length. In the classical limit taken above all 
terms of higher order in 4 (or equivalently in Ti) drop. Of course, the details of these quantum 
contributions depend on the structure of V[S + , S~, S z ], i.e. the ordering of the spin operators. 
Using eqs. (|52"|) , (|53|) expectation values of type (|5l|) can be calculated easily. Here we give a 
list of diagonal elements useful for the derivation of eq. (^0[): 



(55) 
(56) 
(57) 
(58) 
(59) 



S 



2S 



(hSf ( cos 3 + 



2S 2S 2 J 



cos i? sin t? 



(nsy i 



s 



1 



e lLp ( cos 2 sin + — sin 3 (60) 



(S; 0, ip\S + S z S z + S Z S Z S + \S; 0, if) 
= 2{hSfe l * fl ' W 1 



25 



1 



S 



cos 2 sin H sin 



(61) 



(5;0,^|S Z S Z £ Z £ Z |S';0,^) 

= C^ 1 (™* * + (I - w + si) cos2 " sm! * + {w + w) sin! ") (62) 
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(S; tf, tp\S z S z (S+S- + S~S + ) + (S+S- + £-£+)^,S 2 |S; 0, y>) 

= 2(ftS) 4 ((2-1 + ^-^) cos2 ^ sin2 19 + I cos2 * + (| + 4^) sin2 *) (63) 

(5; 0, <p\S z S z {S + S + + SI'S-) + (S + S + + S~S-)S Z S Z \S; 0, y>) 

= 2(hSf cos(2^) ^2 - I + A _ C o S 2 S in 2 + (± - ^ sin 2 0) (64) 

B Coherent states and the classical limit of Spin systems 

Let us consider a system of spins The dynamics are given by a quantum Hamiltonian 

H[{iSi}j e j], which is taken to be an arbitrary polynomial in the spin components. 

We now evaluate the Heisenberg equation of motion for the i-th spin at a chosen time t = to 

in terms of a tensor product of spin-coherent states 

hK*o)> = ® \Si-,0i,<Pi) (65) 



and denote by (•) an expectation value in the state fl65j): 

The right hand side of this equation can be evaluated straightforewardly in terms of the 
variables <pi, while for the left hand side we need information about the time evolution of 
the wave function \ip(t)). 

The time evolution of the spin expectation values is 

(i/>(t)\§iW(t)) = (V(to)| exp(+~H(t - *o))l exp(-~H(t - t ))|^(*o)> 



+ k) t0 + l^j^fm, in, s z ]]) t0 + ■■■ (CD 



The diagonal elements in the expansion (67) are products of expressions of the form ( pl| ) and 
therefore reduce in the limit S — > 00, h — ► 0, ftS = const to the classical values. E.g. for the 
first commutator we have 

hm ~([H,S i ]) t0 = -S i s i x 9{sfgi) (68) 

RS,=S9 ! =const. 

with s*j = (sin #j cos ^ , sin $j sin 99 $ , cos . The right hand side of this equation is nothing but 
the classical Poisson bracket {Sfsi,H} occuring in the well-known Landau-Lifshitz-equation, 
i.e. the equation of motion for a classical spin system. Similar arguments hold for the higher 



commutators in (37), and we end up with 



lim m)\km) = Sfsi + (t- t ){S?s u H} + {t * o) {{Sfs u H},H} + ■■■ (69) 

Sj— *oo 2 

HS i =S'r l =const. 
J 3 
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Thus, we have confirmed that the classical limit of the quantum mechanical time evolution 
( |6?| ) reproduces the motion of a classical spin vector in the case of an arbitrary Hamiltonian. 
By similar considerations one can convince oneself that minimizing the action 

s = £ 2 dtm)\^j t -n\m) (70) 

in the classical limit with respect to the functions $i(t), (fi(t) also leads to the classical equa- 
tions. This should be expected since the variation principle concerning ( |70| ) (with arbitrary 
wave function \ip(t))) is equivalent to the Heisenberg equation of motion. 
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